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The superfluid - Mott insulator phase transition in a Bose-Einstein condensate of neutral atoms with doubly 
degenerate internal ground states in an optical lattice is theoretically investigated. The optical lattice is created 
by two counterpropagating linearly polarized laser beams with the angle 6 between the polarization vectors (lin- 
angle-lin configuration). The phase diagram of the system and the critical values of the parameters are worked 
out. It is shown that the sign of the detuning plays an important role and that there is a strong suppression of 
the Mott transition in the case of blue detuning. Varying the laser intensity and/or the angle 6 one can manip- 
ulate the Mott-insulator to superfluid quantum phase transition as well as prepare the condensate in physically 
distinguishable "ferromagnetic" and "antiferromagnetic" superfluid states. 
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The study of quantum phase transitions (QPT) in optical lat- 
tices has become one of several focusses of current interest in 
the ongoing exploration of the rich physics of Bose-Einstein 
condensates (BEC). In the seminal paper by Jaksch et al. 1 1] 
it was predicted that by increasing the strength of the peri- 
odic optical potential one can confine ultracold atoms at lattice 
sites, which leads to the superfluid-Mott insulator phase tran- 
sition, and this effect has been experimentally observed in a 
one-component BEC of rubidium atoms |2]. Recently coher- 
ent transport of multicomponent BEC in optical lattices has 
been demonstrated 1 3]. In such systems a number of new phe- 
nomena associated with QPT are possible. Theoretical stud- 
ies show that by using different laser schemes one can cre- 
ate coexisting superfluid and Mott phases 1 1, 4|, fragmented 
condensates with topological excitations |5], maximally en- 
tangled atomic states |6], dimer phases |7], and heteronuclear 
molecular BECs 1 8 1 . 

One of the most intriguing features of the multicomponent 

BECs is the structure of their internal levels. The action of 

lasers on multi-level atoms can lead to interference of atomic 

levels. This kind of interference has been used in atomic BECs 
n 

to slow down the light [9], which allows us to excite quantum 
shock waves in a BEC Llfl] . In the present Letter we shall in- 
vestigate the possibility of using interference of atomic levels 
for the manipulation of QPT in a two-component BEC with 
spatially periodic coupling of the degenerate internal ground 
states. It is shown that one can effectively manipulate QPT 
not only by varying the laser intensity, but also the laser po- 
larization, and one can even change the sign of the tunneling 
matrix element, which leads to "feiTomagnetic" and "antifer- 
romagnetic" states in analogy to the spin ordering in magnetic 
systems. 

We consider a two-component BEC of neutral polarizable 
atoms of mass M, possessing an excited electronic state char- 
acterized by the magnetic quantum number m = and two 
Zeeman-degenerate internal ground states with m = ±1 
{Fg = Fe = 1), in a one-dimensional optical lattice. The 
latter is assumed to be created by two counterpropagating 



linearly polarized laser waves of equal amplitudes and fre- 
quencies with the wave number fc^, and the angle 9 between 
the polarization vectors (lin-angle-lin configuration) , and 
with detuning A from the internal atomic transition. In or- 
der to avoid decoherence due to spontaneous emission, |A| 
must be much larger than the spontaneous emission rate. The 
running laser waves form left- and right- polarized standing 
waves with the Rabi frequencies Vt± = fio cos(fciZ ± 9/2). 
Because of the large detuning the excited state can be adiabat- 
ically eliminated. The resulting effective Hamiltonian couples 
the atomic ground states and has the following form; 
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where is a two-component spinor-field operator, (i^ deter- 
mines the periodic potential with the period ir/kL- It is a 2 x 2 
matrix with the components fl^p = ilaflp, a, {3 = +,— . The 
parameter g describes the repulsive interaction of the conden- 
sate atoms. In the one-dimensional case it has the form fl^ 
g = Ah^a/Ma\ (1 — 1.4603a/a^), where a is a symmet- 
ric scattering length and = ^j2fi/Muoi_ is the size of the 
ground state for the harmonic potential with the frequency LOi_ 
confining the BEC in the transverse directions. 

The starting point for the investigation of QPT is the Bose- 
Hubbard model. We assume that the atoms are prepared in 
the lowest Bloch band and the matter-field operator 'J can be 
decomposed in the Wannier basis Ilil3ll 



= ^ exp [iipi) W^(z)ai , 



(2) 



where Wj(z) = W(z — zi) are two-component Wannier 
spinors for the lowest energy band. They are obtained by the 
solution of the eigen-value problem for the Hamiltonian Q 
in the case g = and satisfy the orthonormality condition 
J Wj(z) • Wj(z) dz — 5ij. The indeces i, j label the sites 
of the one-dimensional periodic lattice. The phases are not 
yet defined and their proper choice is discussed in a moment. 
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The Oi and their adjoints are Bose annihilation and creation 
operators attached to the lattice sites. Substituting (|2j into Q 
and taking into account only the hopping between the nearest 
neighbour lattice sites and the atomic interactions at the same 
lattice site, we obtain the well-known Bose-Hubbard Hamil- 
tonian 



Hbh = -J ^ alajexp[i{(p.i - (fii)] 



(3) 



<hj> 



U 



rii ~ 1) ~ ny^uj, 



where /i is a chemical potential. The tunneling matrix element 

= - / Wl+i(z) • (-^^ + ^n') ■ W,(z) dz and the 

atomic interaction parameter U — g J W,| • • W.;^ • 

W; dz can be simultaneously changed by varying the effec- 
tive Rabi frequency and/or the angle 9, but the variation of J 
is much faster Typical dependences of the ratio J/U on the 
dimensionless effective Rabi frequency q — f7§/4ti>/fA, and 
the angle 9 are shown in the left column of Figs.m2j and|3 
(jJh = hk\/2M is a one-photon recoil frequency. The cal- 
culations are performed for ^''Rb (M = 1.45 x 10^^^ kg, 
a = 5.4 nm Q), uj^^2ttx 200 Hz, Al = 27r/fcL = 780 
nm. 

In our model we neglect the antisymmetric scattering 
length, which is responsible for (i) the asymmetry of the in- 
teratomic collisions for the atoms in the same and in different 
ground states and (ii) inelastic collisions of the atoms in the 
states m = — 1 and m = 1, which lead to the creation of the 
atoms in the ground state m = 0. This approximation is well 
justified at least for *^Rb, because the asymmetry leads to a 
correction of the parameter U of the order of only a few per- 
cent in Eq.(|3} and the inelastic collisions become important 
for the times T > 3/n s, where n is the number of atoms per 
lattice site. 

The phases are determined from the requirement of min- 
imal energy of the Hamiltonian (|3jl, which amounts to de- 
manding that J cos{ipi — fi+i) is maximal. In the familiar 
case of a one-component BEC in a periodic lattice, where the 
Wannier spinors reduce to functions, this leads to Lpi = ipj, 
because in such a situation J in Eq.(|3) is real and positive. 
Therefore the superfluid state of a one-component BEC in 
the optical lattice always has a ferromagnetic ordering of the 
phases (pi, which can then all be set to zero. Interestingly, in 
the case of a two-component BEC J can also become neg- 
ative. In this case the choice (pi ~ fi+i does not provide 
a minimal energy and ^pi+i — (pi ± tt has to be chosen in- 
stead. This corresponds to an antiferromagnetic ordering of 
the phases of neighboring lattice sites. Indeed, as one can see 
in Figs.m2j and|3(left column) the parameter J takes in gen- 
eral positive as well as negative values (U is always positive 
for a BEC with repulsive interaction). Therefore, ferromag- 
netic and antiferromagnetic phase ordering occurs in different 
parts of the phase diagram. 

Ferromagnetic and antiferromagnetic phase ordering is pos- 
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FIG. 1: Typical dependences of J/U on the dimensionless effec- 
tive Rabi frequency q — Q,f)/4:L0RA for red detuning (A < 0) and 
corresponding phase diagrams for n — 1. 6 = 50° (a), 64° (b), 
65.7767° (c), 66.5° (d). Horizontal straight lines show the critical 
values (2J/[/)c = ±(3 — 2\/2). Vertical dashed lines on the phase 
diagrams (right column) indicate the boundary between ferromag- 
netic (J > 0) and antiferromagnetic (J < 0) phase ordering. The 
regions of the Mott phase are enclosed by the lobes and loops. 



sible only for the superfluid phase. These two cases are readily 
distinguishable experimentally in the superfluid regime via the 
spatial interference pattern generated by the coherent matter 
waves which one obtains after removing the lattice potential 
1 2]: The interference maxima obtained in the ferromagnetic 
case turn into minima in the antiferromagnetic case and vice 
versa. In the Mott phase, the numbers of particles occupy- 
ing each lattice site are equal and integer, and the phases of 
the corresponding wave functions are completely undefined. 
Therefore, the choice of the ipi remains arbitrary in this case 
and has no observable consequences for the interference pat- 
tern. 

The change of the sign of J is closely related to the form of 
the dispersion relation E{k) for the lowest Bloch band. Using 
the definition of the Wannier spinors the expression for J can 
be rewritten in the form ./ = ^ /g ^ -^'(fc) sin(7rfc/fci) dk. In 
the case of one-component condensate one always has a nor- 
mal dispersion, i.e., E'{k) > for < fc < fc^. In the case 
we are dealing with one can get anomalous dispersion, i.e., 
E'{k) < for < fc < fcL, as wefl. The change of the dis- 
persion type and as a consequence of the sign of J happens at 
the points (q, 9) indicated in Figs.|3land|5l Since the type of 
dispersion is different in the ferromagnetic and antiferromag- 
netic superfluid phases, one can expect principally different 
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FIG. 2: The same as in Fig.Qas functions of S (A < 0). q 
(a), -4.5 (b), -4.80064 (c), -5 (d). 



= -3.9 



properties of the nonlinear atomic matter waves in the two 
regimes IT^ . 

The phase diagram of the Hamiltonian Q is determined 
by the ratio J/U. In the case of one-component BECs it 
is a monotonic function of j^/j. In the two-component case 
we are dealing with, it has quite different properties. Its de- 
pendence on q and 9 is not monotonic, and it can vanish 
or even change its sign at certain finite values of q. There- 
fore, it is reasonable to draw fi — q — 9 diagrams, which are 
shown in the figures, instead of /i — J diagrams. In the mean- 
field approximation and second-order perturbation theory the 
boundary between the superfluid and the Mott regions is given 
by 116] 2\J\/U = {l-n + ^l/U){n^ ^^/U) / [1 + ^l/U), 
wheren — 1 < ji/U < n. The critical values {2\J\/U)^ = 1 + 
2?i— 2 vr?+n, which correspond to {fi/U)^ — V + n— 1, 
define the maximal possible width of the Mott region on the 
phase diagram. Since J/U is not a monotonic function of 
q and 9, there can be several critical values of q and 9 (see 
Figs, n El and|3. At the points [q, 9), where J/U vanishes, 
we have the Mott phase for any values of /i/J7. These points 
define the boundary between the ferromagnetic and antiferro- 
magnetic states. 

The phase diagrams in the plane spanned by 9 and q for 
{li/U)c — \/2 — 1 at n = 1 and negative (red) and positive 
(blue) detuning A are shown in Figs. |3l and [Sj respectively. 
The principal difference between the cases of positive and 
negative A can be understood if we apply the unitary trans- 
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FIG. 3: Phase diagram in the {0,q) plane for {jj./U)c ~ \/2 — 1 
(A < 0). The boundaries between the superfluid and the Mott (n = 
1) phases are shown by solid (calculation for the two-component 
BEG) and dotted (approximation, which does not take into account 
the coupling between the bright and dark states) lines. The dashed 
line corresponds to J = 0. It lies always in the Mott phase, separat- 
ing the regions of the ferromagnetic (SF) and antiferromagnetic (SA) 
superfluid phases. 
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to the Wannier spinors in Eq.(|3}. After the tranformation we 
end up with the bright and dark states iflTil . which are not de- 
generate in contrast to the original ones. The important point 
is that only the bright state is directly coupled to the electro- 
magnetic field and influenced by the potential 



Vb 



A A ^ 



cos 9 cos 2kLz) 



(5) 



We consider first the case 9 = 0, when the transformation 
U does not depend on the position z. In this case the dark 
state does not "feel" any periodic potential. Since in the case 
A > the dark state has a lower energy the atoms stay in this 
state and, therefore, they can not be localized on the lattice 
sites for any laser intensity. In the case A < the situa- 
tion is reversed: The energy of the bright state is lower than 
that of the dark one and only the bright state is populated by 
the atoms. Therefore, increasing the laser intensity, one can 
strongly localize the atoms on the lattice sites in exactly the 
same manner as in the case of the one-component BEC. 

In the case 9 0,U is a position-dependent transformation. 
The atomic center-of-mass motion leads to the gauge potential 



Vn = huJB. 



sinf 



1 + cos 9 cos 2kLZ 



(6) 



acting on the bright and dark atomic states and to the motional 
coupling of the states 1 17|. It seems that the transformation U 
does not lead to any simplifications in the case 9^0. Never- 
theless, it allows one to understand what is going on, assum- 
ing that [Oq/AI ^ Lj/f. In this approximation Vg <C \Vb\, 
and one can neglect the gauge potential for the bright state 
as well as the motional coupling between the bright and dark 



4 




45.6 



45.7 



45.8 45.6 



45.7 



45.8 



FIG. 4: The same as in Figs. Q and |2| for blue detuning (A > 0). 

e = 45° (a), Q = 6 (b). 
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FIG. 5: (a) Phase diagram in the (6, q) plane for {fi/U)c ~ \/2 — 1, 
A > 0. The line J = as well as the two boundaries separating SF 
and SA superfluid phases from the Mott phase are indistinguishable 
on the large-scale plot. The Mott phase is located in the extremely 
narrow region between the superfluid phases. Plot (b), therefore, 
shows the difference Sqi between the two boundaries (solid line) and 
the difference Sq2 between the line J = and the boundary separat- 
ing the ferromagnetic superfluid phase from the Mott phase (dashed 
line). 



states 1 17]. Then the only potentials acting on the bright and 
dark states are given by Eqs. (|5} and (|6}, respectively. On the 
basis of the same argument as in the case ^? = we see that 
the atomic localization in the cases A < and A > is 
determined by the potentials Vb and Vg, respectively. 

Accordingly in the case A < the quantity JIq cos(6')/A 
plays the role of the effective Rabi frequency and defines the 
strength of the periodic potential. Therefore, in the simple ap- 
proximation we consider for the sake of this discussion, the 
critical effective Rabi frequency i^e// — il'^j:j:/2cos9 (dot- 
ted line in Fig. |3, where ^^fj is the critical effective Rabi 
frequency of the one-component BEC and the factor 2 reflects 
the fact that two atomic ground states are involved. Since 



and the transition into the Mott state in the case A > is not 
possible (compare with the exact results in Figs.|3and|5|i. 

As it follows from the results presented in Figs.|3]and|5j this 
simplified description provides a correct physical insight, but 
does not always work and at the laser intensities we are deal- 
ing with one has to take into account the coupling between 
the components, as we have done in the numerical calcula- 
tions leading to the results displayed in Figs. It is also 
necessary to note that the one-component approximation can- 
not describe the change of sign in J and the related transition 
from the ferromagnetic to antiferromagnetic superfluid phase. 

Summarizing, we have studied quantum phase transitions 
of a BEC with two degenerate ground states in an optical 
lattice created by the lin-angle-lin laser configuration. It is 
shown that the periodic coupling of the atomic ground states 
modifies essentially the phase diagram of the system. The 
most surprising feature we encountered are (i) change of sign 
of the tunneling matrix element J under variation of the laser 
intensity and the angle 6, which leads to the ferromagnetic 
and antiferromagnetic phase ordering analogous to the spin 
ordering in magnetic systems, and (ii) suppression of the Mott 
transition in the case A > 0, except for a narrow region in the 
phase diagram (Figs.|4]and|5}. Numerical results presented in 
the figures are for the case of a commensurate filling of the 
lattice with one atom per site. For higher occupation num- 
bers the physics remains the same and there are only minor 
quantitative changes. 
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lOojR |1], the potential Vg appears to be too weak 
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